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Nonperturbative Contributions to the QCD Pressure 
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We summarize the most important arguments why a perturbative description of finite-temperature 
QCD is unlikely to be possible and review various well-established approaches to deal with this 
problem. Then, using a recently proposed method, we investigate nonperturbative contributions 
to the QCD pressure and other observables (like energy, anomaly and bulk viscosity) obtained by 
imposing a functional cutoff at the Gribov horizon. 

Finally, we discuss how such contributions fit into the picture of consecutive effective theories, as 
proposed by Braaten and Nieto, and give an outline of the next steps necessary to improve this type 
of calculation. 
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I. INTRODUCTION 

A. Rise and Fall of the Quark-Gluon Plasma 

One of the most striking properties of QCD is asymptotic 
freedom. For large momentum p, the coupling g(p) is 
small, so quarks and gluons can be treated as if they 
were almost free particles — in particular, they can be 
treated with the sophisticated methods of perturbation 
theory. 

While the situation is obviously different in the regime 
of low energies (which is most relevant for nuclear 
physics), it was natural to expect that a perturbative 
decription could be applied to QCD at sufficiently high 
temperatures. After all, high temperature T implies high 
average particle momentum and thus a small coupling, 
i.e. almost free particles. 

For this scenario, the term "quark-gluon plasma" was 
coined [l[ , and one could expect a phase transition where, 
on a certain curve in the ^-T-space (where fi denotes the 
chemical potential), hadrons melt into such a plasma. 

This phase transition offered a natural solution to a 
problem posed by Hagedorn, 0], who found that due to 
an exponential increase of the number of accessible states, 
the temperature of a hadron could not exceed a certain 
limit T H ~ 160MeV. 

The picture of hadrons melting into a plasma of (al- 
most) free quarks and gluons, however, turned out to 
be too naive, fn principle, this should have been clear 
at least since 1980, when it was shown 0, Q that at 
order g 6 a natural barrier arises for any perturbative de- 
scription. Even earlier than that, the simple fact that 
the infinite-temperature limit of four-dimensional Yang- 
Mills theory is a three-dimensional confining Yang-Mills 
theory could and should have been regarded as a sign 
that any straightforward perturbative approach to high- 
temperature QCD was necessarily doomed. 

It took however more than 20 years until it (slowly) 
began to be accepted that the high-temperature phase of 
QCD has little to do with a conventional plasma. The 
results of the RHIC experiments, [f| , showed clearly that 



also above the phase transition, bound state phenomena 
can not be neglected, and the description as a perfect 
fluid is much more accurate than the one as a weakly 
interacting plasma. 

While the term "quark-gluon plasma" is still widely 
used, one begins to speak (more accurately, though some- 
how using an oxymoron) of a "strongly coupled quark 
gluon plasma" [1, 0] ; or even a "quark-gluon soup" . 

With the experimental results which are — for certain 
observables — an order of magnitude away from the pre- 
dictions for a weakly coupled plasma (see for example 
data on the elliptic flow in Q ) an accurate description of 
the high-temperature phase remains a challenge for the- 
oretical physics. One conclusion, however, seems to be 
clear: In the high-T regime, perturbation theory has to 
be replaced or at least supplemented by nonperturbative 
methods. 



B. Organization of paper 

After the introduction given in HI in |TT] we briefly review 
aspects of finite-temperature QCD. In particular, in III Al 
we discuss which thermodynamic quantities might be in- 
teresting to look at, while in III Bl wc examine the argu- 
ments for a breakdown of perturbation theory. In lH Cl we 
summarize the known perturbative results, which can be 
rederived and extended by effective field theory methods, 
which arc discussed in III Dl In III El we discuss previous 
functional approaches and in III Fl lattice results. In III HI 
we compare the results of these different methods and 
discuss questions of convergence. 

In IIIII we approach the problem with a new "semi- 
perturbative" method [l(| EH which is briefly re- 
viewed in MI Al The physics behind this method is that 
the functional cut-off at the Gribov horizon suppresses 
the infrared components of the gluon field [Hj], so that 
the infrared divergences of finite-temperature field theory 
found by Linde (3[ do not arise [l3|. This method typ- 
ically involves a temperature-dependent renormalization 
scale, an issue we discuss in llll CI In lHI Dl we examine the 
calculational methods used to solve the resulting equa- 
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tions and the expansion used to extract the asymptotic 
form, before we present our results in llll El 

In section IIVI wc discuss these results and how they 
are related to other approaches. In particular, in IIV Al 
we compare different ways to access the nonperturbative 
sector of hot QCD, in IIV Bl we resume the discussion of 
convergence and in IIV CI we present some ideas about 
how to pursue further research. 

In|V]we summarize our results and give a brief outlook. 




FIG. 1: A "ladder diagram" contributing at order g to the 
free energy for I > 3. 



II. HIGH TEMPERATURE QCD 
A. How to Study High Temperatures 

It is useful to rescale thermodynamic quantities with ap- 
propriate powers of the temperature. In particular, the 
free energy per unit volume, 1 the pressure, and the en- 
ergy per unit volume 
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According to up to a perturbative contribution, 
the bulk viscosity C for hot gauge theories is given by the 
logarithmic derivative of the anomaly, 
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where loq denotes a perturbative scale and £y is a per- 
turbative contribution. This formula can be derived from 
the Kubo formula of linear response theory. That the vis- 
cosity is linear in the trace of the energy-momentum ten- 
sor (instead of quadratic) is not surprising in view of 
the Schwinger-Dirac relations, as discussed for example 
in M- 



B. The Perturbative Problem in the Infrared 

Perturbative calculations at finite temperature are dra- 
matically different from those at T = 0. One of the most 
striking differences is that one cannot determine the or- 
der of a graph by simply counting the number of vertices. 



1 In statistical mechanical usage the "free energy" is given by F - 
-wVT. 



Actually a vacuum or propagator graph may be nonana- 
lytic in g 2 . 

While ultraviolet divergences are regulated exactly the 
same way as in the zero-temperature theory, with no ad- 
ditional effort necessary, for T > additional infrared 
divergences appear. They come from Matsubara fre- 
quency n = 0, which has the infrared divergences of 3- 
dimensional Euclidean gauge theory that are even more 
severe than in 4 dimensions. For this reason the Gribov 
horizon, which affects primarily infrared components of 
the gauge field, is more important at finite and high T 
than at T = 0. 

Here, however, another subtlety of thermal field theory 
comes to rescue: Thermal fluctuations give rise to self- 
energy, which, in the static limit p — > corresponds to a 
mass m. At first glance, there are two natural candidates 
for the scale of such a mass, the electric screening mass 
m c \ ~ gT and the magnetic screening mass m mag ~ g 2 T '. 

The mass which is dynamically generated appears in 
the value for ladder dia gram s like the one depicted in 
fig. Q] (see section 8.7 of [l6|). For this type of diagram 
we obtain (ignoring the complicated tensorial structure) 
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for I > 3 . 



If m were independent of g or, like m a \, of order gT, 
we could proceed with perturbation theory without seri- 
ous problems, since an increasing number of loops would 
always correspond to an increasing power of the coupling 
g. It turns out, however, that m is (in the best case) of 
the order of the magnetic screening mass, m mag ~ g 2 T '. 

Thus for any value I > 3 one has contributions of or- 
der g 6 , the perturbative procedure becomes impractica- 
ble unless a suitable resummation technique is available 
— and such a technique has not been found up to now. 



C. Direct Perturbative Approach 

We have seen that the perturbative treatment of the 
QCD free energy runs into fundamental problems at or- 
der g 6 . Still, one can expect that for sufficiently small 
values of g (i.e. for sufficiently high temperatures) the 
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possible perturbativc description (to order g ) still pro- 
vides a good description. 

This is indeed the case (although, as we will see in sec- 
tion [TlIIl only at ridiculously high temperatures). Un- 
fortunatey, even these calculations turn out to be highly 
involved. We summarize here known results, which are 
also collected in [l||, but specialize them to the case of 
pure gauge theory. 

Zcroth order just gives the Stefan-Boltzmann law for 
SU(N) gauge theory, 
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For the second-order contribution, one obtains, [l7|, ]li 
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with Ca denoting the Casimir of the adjoint representa- 
tion, and Ca = N for SU(N). Due to non-analyticity, one 
has a contribution of 0(g 3 ), calculated in [l9j . 
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The g 4 ln<? contribution has been calculated in [2C 
full g 4 term has been obtained in [2l], [22j 
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where 7e denotes the Euler-Mascheroni constant and £ 
the Riemann zeta function. 
At order g 5 , one obtains, [23[ 
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D. Effective Field Theory 

The result of order g 5 is the last one obtained in strict 
perturbation theory. It has been rederived by Braaten 
and Nieto [24| , using an effective field theory method that 
is built on the idea of dimensional reduction 0, [2(| . 

The problem of infrared divergences is adressed by two 
effective theories that are constructed "below" perturba- 
tive QCD. We know that there are three important scales 
present, namely 



2-kT . . . scale of "hard modes" 
gT . . . chromoelectric scale 
g 2 T ... chromomagnetic scale. 

Thus it makes sense to describe each scale in a somewhat 
different way. To do this, two cutoff scales A# and Am 
are introduced, that have to satisfy 
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The region with p > Ae can be reliably described by 
perturbativc QCD, and for this contribution to the free 
energy, called fs, one obtains a power series in g 2 with 
coefficients that can depend on In . 

For Ae > p > Am, with the hard modes integrated 
out, an effective three-dimensional theory, called electro- 
static QCD (EQCD) is introduced, 
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(12) 



where F* = d l A a J -d 3 A a i +g E ! ahc A\A) denotes the mag- 
netostatic field strength tensor and <S£eqcd contains all 
other local (3-dimensionally) gauge-invariant operators 
of dimension three or higher that can be constructed from 
Ai and Aq. The parameters gE, the, are detcmincd 
by matching to perturbative QCD, in particular one has 
m E ~ m e i ~ gT. 

This theory still allows perturbative treatment, making 

2 

use of an expansion in the dimensionless quantities ~ 

g, -£f- etc. This gives for the contribution Jm to the free 
energy a power series in g, with coefficients that depend 
on In j^r and In . The whole series is multiplied by 
the common factor (gT) 3 T. 

The infrared cutoff A M of EQCD is the UV cutoff of 
another theory, magnetostatic QCD (MQCD), 
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with <5£mqcd denoting all gauge-invariant operators of 
dimension 5 or higher. This theory is confining and thus 
truly nonpertubative, but according to [24|, this contri- 
bution to the free energy, called fa, can still be expanded 
in a power series in g, which is multiplied by a general 
factor (g 2 T) 3 . However the value of the coefficient cannot 
be determined perturbatively. 

Since the (well-established) nomenclature may seem 
slightly misleading at first glance, we have tried to give 
a graphical representation in fig. [2j 

MQCD is genuinely nonperturbative, its degrees of 
freedom are (2 + l)-dimensional glueballs. In [26J it was 
suggested to calculate the contributions from this scale 
directly by lattice methods. 

With the effective field theory, it is possible to compute 
the g e ln<7 contribution [27|. The contribution obtained 
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with £if = ff 2 T(l + C(g 2 )) and the constant 

B G = -0.2±0.4 MC ±0.4 SQ , (16) 

where the first error stems from the Monte Carlo simu- 
lation, the second one from the Stochastic Quantization 
procedure employed to obtain the final result. Note that 
Bq = is compatible with this result. 



FIG. 2: The scales of perturbative QCD (pQCD), electro- 
static QCD (EQCD), magnetostatic QCD (MQCD) and the 
different contributions /b, }m and /g to the free energy. The 
coefficients a k and b k are polynomials in logarithms of ratios 
of scales, a k = ft (In £), b k = Q k (ln^, In £-). While the 
coefficients a k and b k can be determined, at least in principle, 
in perturbation theory, this is not possible for c k . 



this way has to be regarded as partly conjectural, since 
the argument inside the logarithm is not clearly defined 
until the full g 6 contribution is known. 

The result thus relies on a supposed structure of can- 
cellation patterns. In addition, it is believed to be reli- 
able only for sufficiently high temperatures (which could, 
however, mean down to T ps 2Tc), since description by a 
three-dimensional theory is valid only for such tempera- 
tures . With these caveats in mind, one obtains for pure 
SU(3) gauge theory 
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with a yet undetermined coefficient qi 0%> for the pure <? 6 
contribution. (See also [H,[29[.) 

This coefficient consists of both perturbative contri- 
butions (from pQCD and EQCD) and nonperturbative 
contributions (from MQCD). It was estimated in (2?| by 
a fit to four-dimensional lattice data for the pressure. 2 

Some of the perturbative contributions of order g 6 are 
known by now (3^, [3l|, but others remain unknown. 
The nonperturbative coefficient has been determined by 
three-dimensional lattice calculations and matching to 
perturbative four-loop calculations in [13] (see also [33[ 
for some cases with N ^ 3) and [34|. One obtains 
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2 The problem with such a procedure is that in the regime where 
lattice data is available, the contributions of higher order may 
also be large. 



E. Functional Approaches 

Due to the limitations of perturbation theory, nonpertur- 
bative methods definitely deserve a closer look — more- 
over the effective field theory approach also relies on the 
ability to calculate certain quantities nonperturbatively. 

As in the zero-temperature case (HJ, also for finite 
temperature, fundamental aspects of Yang-Mills theory 
and QCD are accessible to functional methods based on 
Dyson-Schwinger equations (DSEs) , [H, [3?| . 

For certain asymptotic situations (deep ultraviolet, 
deep infrared, infinite temperature limit) several analytic 
results can be obtained; but in general numerical studies 
of truncated DSE systems are necessary. 

In addition to the standard truncations, finite tem- 
perature calculations require also some treatment of the 
Matsubara series. Usually it is replaced by a finite sum, 
even though this means that the limit of four-dimensional 
zero-temperature theory is now technically hard to ac- 
cess. 

At the present level, these restrictions make it diffi- 
cult to obtain precise quantitative results. Nevertheless 
there is reasonable confidence about the qualitative pic- 
ture that arises from these studies. Both from infrared 
exponents and from numerical results one sees that the 
soft modes are not significantly affected by the presence 
of hard modes, thus the confining property of the theory 
cannot be expected to be lost in the high-temperature 
phase. 

Consequently, while the over-screening (which would 
attribute an infinite amount of energy to free color 
charges) of chromoelectric gluons is reduced to screening 
(as it is the case for electric charges in a conventional 
plasma), chromomagnetic gluons remain over-screened 
and thus confined, which renders any description of such 
gluons as almost free (quasi-)particles meaningless. 

While the functional method yields considerable in- 
sight into propagators and related quantities, unfortu- 
nately the pressure (and quantities derived from it) are, 
up to now, difficult to access in this approach. Neverthe- 
less the results obtained so far by functional methods pro- 
vide additional evidence for the picture of bound states 
playing an important role even at very high temperature 
and part of the gluon spectrum (the chromomagnetic sec- 
tor) being confined at any temperature. 
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FIG. 3: Rescaled pressure of SU(3) lattice gauge theory, 
from [38j], where T c is the transition temperature, and N T = 
4,6,8. 



F. Lattice Gauge Theory 

Lattice Gauge Theory is generally considered the most 
rigorous approach to nonperturbative QCD, and so it is 
natural to also study thermodynamics on the lattice. 

The drawbacks of the method, however, are known as 
well: To reliably approach the thermodynamic and the 
continuum limits, extrapolations which require calcula- 
tions with various different lattice sizes are necessary. 
The inclusion of fermions is expensive, especially if good 
chiral properties are required. 

Despite these drawbacks, lattice data is (apart from 
possible experiments) the thing one usually compares any 
other calculation to. For pure SU(3) gauge theory, the 
problem of determining the equation of state is regarded 
as solved since the publication of [38|, where results are 
further confirmed in [39| . (Note, however, that there 
remain certain doubts about the accuracy of the infinite- 
volume limit, see [40|.) The results for the pressure and 
the anomaly are displayed in figs. [3] and [4] 



G. Other approaches 

Pisarski [4l[ has observed from lattice data [38[ that (e — 
3p) /T 4 x T 2 is approximately constant in a broad range 
above T c . From this and physical reasoning he obtained 
the formula 

Ppure glue(T) « / P crt(T 4 - T/T 2 ). (17) 

It is notable that no T 3 correction is apparent. 
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FIG. 4: Rescaled anomaly of SU(3) lattice gauge theory, 
from [H 

An active area of research is the AdS/CFT duality (4^ | 
and AdS/QCD duality [43|, including in particular du- 
ality at finite temperature (for a pedagogical introduc- 
tion see [HI)- In this connection it is interesting to note 
that a formula similar to Pisarski's has recently been ob- 
tained from this duality and the (truncated) entropy den- 
sity of the horizon of a deformed Euclidean AdSs black 
hole, 0. 



H. Comparison of Results 

Knowing that perturbation theory is limited to some 
fixed order in g, we can still estimate how good the pos- 
sible perturbative description actually is. Ways to judge 
this are to check whether contributions from higher or- 
ders are small compared to those from lower orders or 
to compare perturbative expressions to results of lattice 
calculations. 

Unfortunately, both methods suggest that the conver- 
gence is extremely poor for temperatures of the order of 
several T c , where T c is the transition temperature, and 
to obtain good conver gen ce one has to look at least at 
the electroweak scale, [23l |27| . A plot of the results of 
optimized perturbation theory are given in fig. [5] 

It has been conjectured, [13], that the results of order 
<7 6 are not significantly changed by higher orders (since 
one can hope to have obtained at order g e the main con- 
tribution from each scale; perhaps also due to the fact 
that originally large terms of higher orders cancel against 
each other). To the knowledge of the authors, however, 
there is no strong evidence to support this conjecture. 
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FIG. 5: The convergence of an (optimized) perturbation se- 
ries for "long-distance contribution" to the pressure, from [4(|. 
The order of the expansion is characterized by the dimension- 
less parameter y ~ 9 J , where g\ denotes the gauge coupling 
in the effective three-dimensional theory. The perturbative 

(2 \ 3 
T ) ' 

inclusion of lower powers of y corresponds to higher orders in 
perturbation theory. 



From the existing data one cannot even exclude the 
unsettling possibiliy that for "physical" temperatures the 
perturbation series already begins to diverge at some or- 
der n < 6. This would mean that contributions from 
higher orders are of comparable magnitude to those of 
low order and no systematic cancellations occur. If this 
were indeed the case, we could not expect to have any 
reliable perturbative description for temperatures which 
are accessible in current experiments. 

One should mention that there is an additional ambiguity 
in the perturbative results. All terms beyond the Stcfan- 
Boltzmann contribution contain some power of the run- 
ning coupling g. Thus, for all practical calculations there 
is some dependence on the scale /i, at which g(fj.) is eval- 
uated. Traditionally, one chooses /u(T) = 2irT in the 
high-temperature regime, but there are alternative ap- 
proaches, for example a ppl ication of the principle of min- 
imal sensitivity, [48l I491T50I] . We will discuss this question 
in more detail in section Ull CI 



III. THE SEMIPERTURBATIVE APPROACH 
A. Equation of State from a Local Action 

An alternative approach that combines nonperturbative 
elements with perturbative expansions has been devel- 



oped in 0, [13, HH, which we now briefly describe. 

The basic physical idea is that the infrared divergences 
of finite-temperature perturbation theory do not arise 
when the domain of functional integration is cut-off at 
the Gribov horizon. The cut-off will be done in Coulomb 
gauge which is well adapted to finite-temperature calcu- 
lations. Indeed both the gauge condition, diAi(x, t) = 0, 
and the cut-off at the Gribov horizon are applied to 3- 
dimcnsional configurations on each time slice t, and are 
entirely independent of the temporal extent of the lattice 
< t < P, where f3 = 1/kT. 

The functional cut-off at the Gribov horizon is effected 
at first by adding a non-local term Snl (A) to the action 
[5ll . I52I ]. The non-local term then gets replaced by a 
local, renormalizable term Sl in the action by means of 
an integration over a multiplet of auxiliary fcrmi and bosc 
ghost pairs , 

exp[— Snl{A)\ = / dipdipdcodui exp[— <Sx(y, <p, w, w)]. 



(18) 

The BRST symmetry is explicitly broken by this term, an 
effect which, alternativel y, m ay be interpreted as spon- 
taneous BRST breaking [53]. Although the breaking of 
BRST invariance precludes the definition of obscrvables 
as elements of the cohomology of the BRST-opcrator, the 
eqivalence to the canonical formulation has been estab- 
lished @ , thereby ensuring the physical foundation of the 
approach, including unitarity. Here the physicality of the 
Coulomb gauge plays an essential role. 

The new term in the action depends on a mass param- 
eter m which appears in the Lagrangian density 
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Di(<p - <p)i + g{DiC x uji) 



The adjoint part of the Bose ghost (ip — (p)i mixes with 
the gauge field A; through the term Di(ip — (p)i = 
[pi + gAix)(ip — (p)i. At tree level one obtains a gluon 
propagator, 
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that satisfies the Gribov dispersion relation 



(20) 



(21) 



The functional cut-off at the Gribov horizon imposes 
the condition that the free energy W or quantum effective 
action T be stationary with respect to m, 



dW 
dm 



dm 



0. 



(22) 



This "horizon condition" has the form of a non- 
perturbative gap equation that determines the Gribov 
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mass to = m(T, Aqcd), and thereby provides a new 
vacuum, around which a perturbativc expansion is again 
possible. 

The most powerful non-perturbative methods available 
are called for to solve this system. However in the present 
work we shall modestly investigate a semi-perturbativc 
method Q , in which one calculates all quantities pertur- 
batively in g, including T, taking to to be a quantity 
of order g°, and then one substitutes for to the non- 
perturbative solution to the gap equation (|22|) . We shall 
find that this method can be a good approximation only 
at extremely high energies. Nevertheless as a matter of 
principle, it is a significant success that for thermody- 
namic observables this procedure gives finite results pre- 
cisely at the order, g e at which ordinary perturbation 
theory diverges. 



B. The Gap Equation 

In lowest non-trivial order in the semi-perturbative 
method Q, the gap equation (j2"2")) reads after separation 
in an TO*-dependent and a T-dependcnt part 



f(m*)=y(T) 

OC 

/(to*) := iln — + f 
J K ' 2 to* J 
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where to* = 



__i ln ejm 

Ng 2 (fi) 4 2T 2 



(23) 
(24) 

(25) 



m/T is the rescaled Gribov mass and 
u(x) = \x 2 + \ (26) 



is the reduced dispersion relation. An important source 
of ambiguity, shared with other (semi)perturbative ap- 
proaches, is the choice of the temperature-dependent 
scale n(T) at which the coupling g is evaluated. 



pure gauge theory). Actually, as long as we are in the 
confining region (i.e. below T = T c ), the mass of bound 
states will always be "more important" than the thermal 
energy. 

These restrictions, together with some conditions of 
"naturalness" , can be summarized by demanding that 
the renormalization scale /i(T) should fulfill: 

(I) fj,(T) w fi = const for T < T c , 

(II) ^(T)w27rT for T > T c , 

(III) continuous, 

(IV) monotonically rising for all T, 

(V) convex for all T 

Conditions (I) and (II) might be replaced or supple- 
mented by the asymptotic conditions 



(I') nm( M (T)- / io)=0 ) 



limM (r)- Mo = 

T^Q T 



(IF) lim r (n(T) - 2irT) = for all n e N Q 

T — ►co 

(27) 

Due to the phase transition at T = T c ks Aqcd, a simple 
choice is 



Mo,i(T) : 
Mo,h(T) 



2ttA 
= 2ttT 



MS 



for T < Ajjg 
forT>A^. 



(28) 



This choice is supported by the fact that 27rA-^jg is in 
the order of magnitude of glueball masses. Another rea- 
sonable choice is 



fx(T) = 2nT 



27rA MS C 



(29) 



This form, however, is less favorable for numerical rea- 
sons, thus we have exclusively used (|28|) in the numerical 
studies performed in section fill Dl Both forms are plot- 
ted in fig. 6. 

Another ansatz, used especially in functional calcula- 
tions, [37j ; is a 't Hooft-like scaling 



C. Choice of the Renormalization Scale 

We consider the coupling g 2 (//) at some renormalization 
scale fJ-(T). For a certain temperature T, the optimal 
renormalization scale should be chosen equal to the scale 
that governs the behaviour of the system. For field theory 
at high temperatures, this scale is expected to be equal 
to the lowest Matsubara frequency, i.e. 27rT; for small T 
it should be constant. 

Since we are considering a confining theory with a mass 
gap, for low temperatures the optimal renormalization 
scale is not expected to go to to zero. For a system at 
very low (even zero) temperature, the most characteristic 
scale is not the very small average kinetic energy, but 
instead the mass of the lightest physical object, which is 
some bound state (a hadron in full QCD, a glueball in 



ff 2 (/i(T))T = const, 



(30) 



which, at one-loop level, corresponds to an exponential 
growth /i(T) = fio c aT with some positive constant a. 
This choice provides a smooth infinite-temperature limit, 
but does not respect condition (II), and has not been used 
in the current article. 



D. Calculational Methods 



The gap equation (|23[) is an implicit equation for to*(T), 
which, in constrast to "genuine" integral equations, can 
be solved independently for each temperature T. Our re- 
sults have been obtained in Mathcmatica by combining a 
numerical equation solver with adaptive Gaufi-Lcgendrc 
integration. 



8 



20 
15 
10 
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■ 2/1 

■ 2ttT 



2/1 
2/IT 

fj/T 



FIG. 6: The choices (|28|l and (I29|) for the renormalization 
scale n(T), displayed together with the asymptotics /^o,i = 
27tAqcd, i-to,h = 2ttT. In addition we show and fj,' — 4^. 



The derivatives necessary to obtain anomaly and bulk 
viscosity (see equations (|45|) and l[50])) can be done ei- 
ther numerically or analytically. The second way unfor- 
tunately involves additional integrals, which can again 
only be evaluated numerically. (See appendix for de- 
tails.) 

While both methods are potentially susceptible to nu- 
merical problems, they are of very different nature. Ac- 
tually, the results of both methods agree remarkably well, 
inspiring confidence in the stability of the result. 

All calculations directly involving T have been per- 
formed on logarithmic temperature scale. This allows di- 
rect implementation of logarithmic derivatives, reduces 
numerical errors as compared to calculations on a linear 
scale and enables one to reach significantly higher tem- 
peratures. 

For all quantities under consideration, we could obtain 
asymptotic expressions by expansion in the coupling g 2 . 
In general, we use, [30j j . 



2 ^ op 26 ln " 



? 2 (M) 
dg 2 



b ° lnf21n- /i 



MS 



(31) 



a/j, 



(47T) 



with the coefficients 
00 = -2 (iir) 2 b 



2 , 2-loop 00 4/ \ , 01 

0*) + 
22 C A + ST f 



g\iA (32) 



N 



pure SU(3) 



-22. 



(33) 



1 = -2(^)H 1 = 

pureJU(3) _^ 



-68C A + 4QCa7> + 24C f T f 
N 



(34) 



and the group-theoretical factors Ca = N and Cf = 
N 2 — 1. Ty is equal to half the number of quark flavors 



(a) 



Am* 




10 20 30 40 50 60 



ln- 



FIG. 7: The rescaled Gribov mass m* = ^r: (a) solid - nu- 
merical solution, dotted - asymptotic expression from (I35|l ; 
(b) relative deviation Am* ol = (m^ uln — m aay )/m asy 



and thus vanishes in pure gauge. While the results in 
subsection IIII El are given for the one-loop form (easily 
obtained by setting 0\ = bi = in (f3Tj) and l[32|)). there 
are only minor changes when switching to the two-loop 
form. 



E. Results 

We now summarize the results obtained by numerically 
solving the gap equation and the corresponding asymp- 
totic expressions. 



1. Gribov Mass 

Solving the gap equation yields the Gribov mass m(T). 
An expansion gives to leading order in g 2 



m*{T) 



N 



2 3 / 2 3 7T 



5 2 (M) 



(35) 



The numerical result and this asymptotic form are dis- 
played in fig. [7] The agreement is excellent down to the 
phase transition (below which the formalism is probably 
not applicable anyhow) , thus higher-order corrections to 
the Gribov mass are small. 
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FIG. 8: The rescaled reduced free energy w r — u> r ,sB: (a) 
solid - numerical solution, dotted - asymptotic expres- 
sion from (|38p ; (b) relative deviation Aw ItTC \ — (w Itnum — 

U) r ,asy)/(u'r,asy — Wit.Sb) 



2. Free Energy and Pressure 



For the pressure p and the free energy w we obtain 



J>4 



T3 



(N 2 1) 



1 



2N g 2 {^) 3tt 2 T 4 



K{m) 



K(m) 



dk k 4 



E(k) e^« - 1 



(36) 
(37) 



with k = \k\. An expansion for K(m) is not completely 
straightforward due to a nonanalyticity in m 4 , but, as 
shown in it can be performed and yields the asymp- 
totic expression 



(iV 2 - 1) 



5 (Ng 2 



45 



18 V 4tt 2 



(38) 



The full solution and the asymptotic form are given in 
hg. [51 where we have subtracted the Stefan-Boltzmann 
part, denoted by wsb- In contrast to the case of m, 
higher-order corrections are obviously not small for w 
since agreement between the full (numerical) and the 
asymptotic result is not good below T ss 10 6 A A j s . 

It is instructive to see that K{m) can also be evalu- 
ated by using an intermediate cutoff. While more cum- 
bersome, this method allows us to identify contributions 



from different scales and thus gives some idea how to re- 
late this result to the one obtained by effective theory 
approaches (see section Hi Dp . 

To do this, we introduce a cutoff A with m « A « T, 
which separates contributions from the scale m ~ g 2 T 
and from the scale 2-kT. Doing so, we obtain 



K{m) 



E(k) e^W - 1 J E{k) e^( fe ) - 1 : 

A 



K 2 



where 



(39) 



A/r 



dx 



A/m 



u{x) (e^ 3 "" 1 ^) - 1) 



A/m 



m T 



m 6 T 




The first integral is trivial; for the second we can replace 
the upper limit A/m by oo and apply residue calculus to 
obtain 



T 



A d - 



V2 



rrfT . 



(41) 



For K 2 we obtain 
K 2 r 



dk k 4 4 f 

k eP h -l V ev-1 



T 4, 



A/T 
A/T 



dy 



cv - 1 



dy 



ey - 1 



(42) 



The first integral is the well-known Planck integral. In 
the second one we can again expand the exponential, 
since j/ < | < 1, and obtain 



A/T 



K 2 « T 



' 4 71 — — A 3 



dyy 2 )=T 4 



15 



(43) 



This gives 



15 y/2 



(44) 



The cutoff-dependent parts in K\ and K 2 precisely can- 
cel, leaving a clear separation of the Stefan-Boltzmann 
contribution from k ~ T and the contribution from the 
scale k ~ m ~ g 2 T. 
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3. Energy and Anomaly 

From the free energy or the pressure, we can calculate 
the rescaled anomaly via 

d 



= e - 3p ^ T _d_P_ 
T 4 dT T 4 



P 



d(ln£)T 4 



(45) 



(with some arbitrary scale A) since from ([T]), we have 



T —JL = T— — 

dT T 4 ~ dT T 3 

_ 1 dT 



_ dT 

3wT 



- 3 



T 3 

3 P 



T 4 " T 4 • ( 46 ) 

It is also obvious that the energy can be directly obtained 
from the anomaly by using the relation e = 3p + A. Thus 
we do not show separate graphs for e. From (|45l) it is 
also clear that a// deviations from the Stefan-Boltzmann 

4 

pressure p rj sB = fg- are encoded in the anomaly, since 
integration gives 



Pi(T) =Pr,SB 



MT<) 



dT' 



From ([38]) . (|45|). and ((32]) we obtain 

TV 3 



A r - -(iV 2 - 1) 



3456 7T 4 



(47) 



(48) 



for the asymptotic expansion. The numerical result and 
the asymptotic form are shown in fig. [9] Again higher- 
order corrections are large except for extremely high tem- 
peratures. 



4- Bulk Viscosity 

In formula (UJ) there is one ambiguity, the choice of the 
scale tt>o. According to [54| a reasonable range of val- 
ues is u>o = (0.5 -j- 1.5) GeV. Neglecting the perturbative 
contribution from Ey , we obtain 



J_T 5 —- 
9wq dT 



3 P 



rp4 



The rescaled bulk viscosity is given by 
9 dl 



(49) 



(50) 



In the asymptotic expression, correction terms origi- 
nating from — become quite complicated. Since they 
are relatively unimportant for reasonable choice of /i(T) 
(and even vanish identically for the simple form (|28p ) we 
only give the simplified expression, where we have set 

Z d£. — 1 

H dT ~ x ' 



Cr 



1 N 3 (N 2 
9wo 



1) 



3456 7T 4 
2<rV) P(g 2 )- 



? 4 (m) 



dg 2 



Pig 2 )- (5i) 




FIG. 9: The rescaled anomaly A r : (a) solid - numerical so- 
lution, dotted - asymptotic expression from (148 ) ; (b) relative 



deviation AA r , r ci = (A T 



Ar 5 asy ) / j4r,as\ 



The full expression is derived in appendix [B] Graphs for 
the numerical solution and the asymptotic expression are 
shown in fig. [10] for the choice ujq = 5Aj^g. 

The behaviour close to T = A A j s is strongly influ- 
enced by the choice of /i(T). Apart from that however 
the viscosity £ r rises significantly when the temperature 
approaches the critical temperature from above, in agree- 
ment with 1541. 



IV. DISCUSSION 

A. Access to the nonperturbative sector 

Various results make clear that finite-temperature QCD 
contains in principle a perturbatively accessible sector, 
which, starting at order (g 2 T) 3 , interacts with a gen- 
uine nonperturbative sector. At least formally an ex- 
pansion in powers of the coupling g is possible also for 
non-pcrturbativc contributions. 

According to Gribov's confinement scenario, [H, [5f| 
|56| , the vicinity of the Gribov horizon dominates the non- 
perturbative aspects of the theory. So correctly taking 
into account this region should give access to the non- 
perturbative sector of the theory at high temperatures 
also. Indeed, the cutoff at the Gribov horizon employed 
in this article gives a finite nonperturbative contribution 
to the free energy at order g 6 , where the nonperturbative 
sector of the theory begins to spoil direct perturbative 
approaches. 
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FIG. 10: The rescaled bulk viscosity: (a) solid - numerical so- 
lution (r,num, dotted - asymptotic expression Cr.asy from 
(b) relative deviation AC r , rc i = (Cr.num — Cr,asy)/Cr,a Sy 



The nonperturbative sector (described by MQCD in 
the picture of section IIID|) is also accessible to lattice 
calculations. Comparison of our analytic result (f55|) with 
the lattice expressions (fT5|) and ([T6]) gives 



(6) 

np, analyt 



(N 



1)N 3 f 

— —^gT 3 

10 368?r 4 * 



w 



(6) 

np, lattice 



(iV 2 -l)iV 3 
1 280 7T 4 



(l±4)g 6 T 3 



(52) 
(53) 



These results are compatible, though the errors of the 
lattice calculations are too large at the moment to allow 
a definite statement about the quality of agreement. 




FIG. 11: The rescaled free energy in the low-temperature 
region (solid - numerical solution, dotted - asymptotic ex- 
pression from H48j l) 
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— . ln-^- 
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FIG. 12: The rescaled anomly e ~4 P in the low-temperature 
region (solid - numerical solution, dotted - asymptotic ex- 
pression from Q48| l1 



poor or that there is even no convergence at all. The lat- 
ter would not be completely unexpected. It is well known 
that in quantum field theory series obtained by expan- 
sion in the coupling are rarely convergent, but at best 
asymptotic (and usually even this cannot be proven). 



B. Convergence of the series? 



As already mentioned in section III HI the convergence 
of perturbation series is extremely poor for temperatures 
O(GeV) or below. As discussed in [4{j], this can be traced 
back to the poor convergence of contributions from the 
EQCD sector, which begin to contribute at order g 3 . 

A similar behaviour seems to be true for the contri- 
bution from MQCD. While a formal expansion in g is 
possible (and for very high temperatures T > 10 10 GeV 
the agreement is reasonably good) , the expansion has lit- 
tle to do with the full result for low temperatures. From 
the low-temperature graphs displayed in figs. [11] to [T3] 
it is likely that higher-order corrections cannot be small 
compared to the leading term. 

This suggests that cither the convergence is extremely 
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FIG. 13: The rescaled bulk viscosity in the low-temperature 
region (solid - numerical solution, dotted - asymptotic ex- 
pression from H48j l) 
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Assuming that expansion in g of the QCD free energy 
yields a divergent asymptotic series would give the fol- 
lowing scenario: For each temperature T one can expect 
that an "optimal order" n exists beyond which the series 
leaves the "path of apparent convergence" . For low tem- 
peratures and thus large couplings this order may be so 
small that no partial sum of the perturbation series can 
serve as a satisfactory approximation. 



C. Further Steps 

Our studies leave open several questions that are worth 
further investigation. While the analytic result (|52|) is 
compatible with the lattice result (|53j) . a reduction of 
the errors of the latter would be highly desirable. (Un- 
fortunately the errors would have to be reduced about 
one order of magnitude, a task which would require con- 
siderable computational resources.) 

On the analytic side, higher-order calculations in the 
semi-perturbative formalism could help to further clear- 
ify the connection of our approach to the sequence of 
theories discussed in subsection III Dl They could also 
help to reveal if (j52|) is indeed the full contribution to 
order g 6 from the magnetostatic sector or if there are 
additional contributions from (formally) higher orders as 
well, which are not present in the lowest-order approxi- 
mation to the gap equation (|23l) . 

Such calculations could also shed some more light on 
the question of (apparent) convergence, as just discussed 
in subsection IIVBI Of course also determination of the 
full g e contribution to the free energy (conceptually pos- 
sible in the framework of effective theories) would be very 
helpful for further statements about convergence issues. 

Our results indicate that the semi-perturbative method 
of calculation is reliable only at extremely high T. A 
more advanced approach to the nonperturbative sector 
could involve solving the Dyson-Schwinger equation for 
the system with local action and auxiliary fermi and bose 
ghost pairs Q or, alternatively, studying bound-state 
equations for glueballs in MQCD. Those objects, which 
determine higher-order contributions from this sector are 
closely related (though not strictly identical) to chromo- 
magnctic glueballs in the four-dimensional theory. 



method provides a framework for purely analytical stud- 
ies beyond the limits of thermal perturbation theory, 
without the necessity of lattice calculations (neither in 
the full nor in an effective theory.) 

While for each quantity we were able to obtain the 
leading coefficient of the expansion in the coupling g, we 
also noticed that higher-order corrections have to be large 
in order to accomodate the numerical results. This is one 
more sign that, in thermal QCD, expansion in the cou- 
pling cannot be expected to give reliable approximations 
except for extremely high temperatures. 

While results are promising, further investigations are 
necessary to clearify the limits of this method, illuminate 
further the relation to other approaches and perhaps ex- 
tract general information about the quality of series ex- 
pansions in thermal QCD. 
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APPENDIX A: THE ANOMALY FROM AN 
ANALYTIC DERIVATIVE 

We now show how the anomaly can be obtained without 
the necessity for numerical differentiation. In the follow- 
ing, we always understand g 2 = g 2 (n) and [i = n(T). 

One-loop expansion of the gap equation gives for the 
rescaled pressure 



V. SUMMARY AND OUTLOOK 

We first presented a short synopsis on current problems in 
thermal QCD, where we also reviewed several methods 
to approach them, including consecutive effective field 
theories, as proposed in [24j]. 

Then, using the local action proposed in Q, we ob- 
tained nonperturbative contributions to several thermo- 
dynamic obscrvablcs, including free energy, anomaly and 
bulk viscosity. Being directly related to other non- 
perturbative approaches (like magnetostatic QCD), this 



p _ w 



(N 2 - 1) 



3 TO* 



2N g 2 {n) 3tt 2 



K(m* 



with K{m*) — TO* 4 /(m*) and 



/(to*) 



da; x — 1 



u{x) e m * u <» - 1 ' 



(Al) 



(A2) 



using the reduced dispersion relation l|26p. From this, we 
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obtain 



Diffcrcntion of the right-hand side yields 



A X = T — 

dTT 4 

= {N 2 - 1) 
W 2 -l 



3 dm* 4 1 dK(m* 
■ T — + — - T ■ y 



2N dT g 2 3tt 2 



dT 



2N 



m* 4 do 2 4 m* 3 dm* 

T a T 

g* dT + g 2 



dT 



N 2 - I dK{m*) 



3tt 2 



dT 



(A3) 



The derivative of g 2 with respect to T be easily calculated 
using the /3-function, 



d(g 2 ) _ T d(g 2 ) _ T fl , 



dfi [i d/i 



Kg 2 )- 



(A4) 



dy(T) 

dT 



d 

dT 



3^ 2 



Ng 2 (v) 
3^ 2 d<? 2 (/i) 



1 e^T) 
4 2T 2 



Ng*Qi(T)) dT 
3tt 2 1 

1 T 1 / d ^ 
1 J [ dT 



AI ln 

dT2 

dff 2 (/i) 
dT 



ii 

T 2 



2 yu(T) \ T 
1 ( 3n 2 (3(g 2 ) d f i 



At(T) |_ AT g 4 (fi) dT 2\dT T 



e/i(T) 
2 T 



1 / d/t /i 



(A13) 



From this, we obtain 



dg 2 {v) _ T d 5 2 (/i) d/i(T) _ fl/ a „ T d/i(T) 



dT /x ' d/i dT 



For the last term in (IA3I) we have 



A* dT 



(A5) 



T dK(m*) _ T dm* dK{m*) 



dT 
8K(m*) 
dm* 



dT dm* 



(A6) 



4m* 3 /(m*) -m* 4 T(m*) (A7) 



L(m*) = / dx 



(x A - 1) e" 



i(x) 



(A8) 



We obtain T by differentiating the gap equation 
/ (m*) = y(T) with respect to T, where 



dx 



1 



f (m*) = — In h , 

y ' 2 m* J u(x) e m »W - 1 
o 



3tt 2 _l ln C^ 2 (T) 



7Vg 2 (Ai) 4 2T 2 



(A9) 
(A10) 



Since the left-hand side depends on T only implicitly via 
m*, we have 



d f, *\ d / dm* 

dr /(m } = d^ 



2 m 



J(m* 



dx 



^m M u(x) 



dT ' 
(All) 

(A12) 



The second term inside the curly brackets is a measure 
for the deviation from the asymptotic behaviour; in the 
asymptotic regime with /i(T) — 2nT, the above expres- 
sion simplifies to 



dy(T) 



dT 



asympt 



3tt 2 1 /?(g 2 (2^T)) 
' N T g 4 (2irT) 



(A14) 



Collecting our results, we obtain 



T ■ 



dm* 
~dT~ 



T 



dT 



1 



Ira 
T 



M^) ^ F + J{m*) 
3tt 2 /3(g 2 ) d/i 1 
JV 9 4 dT 2 



9/i /i 
dT~ T 



(A15) 



and the rescaled anomaly is given by 



A r = 3 



N 2 -l 



2N 

N 2 - 



2 T d/i 4 1 



— T ^(ff ) - 377; 



■T- 



dm* 



3tt 2 



- (4m ,3 J(m')-m ,4 L(m*)) T 



dT 

dm' 



dT 
(A16) 
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APPENDIX B: ASYMPTOTIC BULK VISCOSITY 



Keeping the full dependence on /x(T), we obtain for the 
bulk viscosity (|49|) . again making use of (|A4|) . 



f= 1 T 5 d e-3p 
^ 9w dT T 4 

1 iV 3 (TV 2 - 1) 
9^ 3456 7T 4 

1 iV 3 (N 2 - 1) 
9^ 3456 7T 4 



d_ 

dT' 



rpQ 



2g 2 (v) /3(ff 2 )+3 4 (M) 



dff 2 



(j, dT 



d(g 2 ) d/x T d/x 



g\»)P(g 2 ) 



1 dxt 
/x dT 



T /dM 

2 I dT 



/' 



d/x dT /x dT 

T <i 2 /j. 
^dT 2 



1 N 3 (N 2 - 1) 
9^ 3456 7T 4 



2f3(g 2 )+g 2 ^) 



dP(g 2 
dg 2 



T dxi 
/x dT 



T d/x 
/x dT 



T d/x 
/x dT 

2 



i d /x 
x7 dT 2 



(Bl) 
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